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^T^The  design  of  a combined  passive  and  active  sonar  system  in  the 
form  of  a programmable  special  purpose  computer  is  feasible  with  current 
technology.  In  fact,  the  flexibility  offered  by  such  a system  brings 
to  the  foreground  an  important  sonar  problem:  What  algorithms  are  most 


effective  in  accomplishing  a given  task  or  mission?  Since  many  of  the 
algorithms  are  non-linear  and  may  not  be  directly  analyzed  in  terms  of 
relevant  performance  parameters,  a Monte  Carlo  simulation  has  been 
developed.  The  essence  of  the  technique  is  to  develop  a representative 
stochastic  sequence  and  carefully  observe  the  algorithm  (or  algorithms ) 
response  under  controlled  conditions. 


The  Soft1 


The  aim  in  developing  the  computer  simulation  programs  was  to 
build  in  enough  flexibility  to  provide  for  the  simulation  of  several 
different  types  of  signal  processing  hardware  and  to  allow  either 
synthesized  signals  and  noise,  (Monte  Carlo  Techniques)  or  taped  sea 
data  as  input.  By  allowing  for  modification  to  key  programs,  direct 
system  comparisons  can  be  made  with  identical  stochastic  samples.  The 
reproduc ibility  and  exact  control  of  synthsized  signal  and  noise 
characteristics  are  the  most  important  features. 


Detailed  in  thi6  report  are  six  specific  computer  programs  which 
are  the  result  of  a multibeam  sonar  study  in  progress  by  Code  2242. 
Hie  six  programs  are  shown  in  Figure  1 and  will  be  detailed  in  that 
order.  The  first  three  programs  are  used  to  synthesize  a realistic 
acoustic  field,  and  the  last  three  are  used  to  simulate  the  hardware 
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and  extract  certain  statistics. 

Synthesis  of  Independent  Filtered  Random  Noise 

The  purpose  of  the  first  program  was  to  produce  several  channels 
of  independent  filtered  random  noise  with  a specified  spectral  shape. 

The  noise  is  produced  by  initially  generating  a series  of  independent 
and  normally  distributed  random  numbers  and  then  "providing  time 
correlation  or  dependence"  with  digital  filtering  techniques  to  produce 
the  desired  spectral  characteristics  for  each  channel. 

The  techniques  used  to  generate  independent  but  repeatable  random 
numbers  by  a digital  computer  are  detailed  in  an  excellent  survey  paper 
by  Chambers1.  Besides  describing  the  history  and  application  of  random 
number  generation  by  computers,  a survey  was  made  of  current  algorithms 
and  their  shortcomings.  Chambers  conclusion  was  that  perhaps  the  best 
algorithm  is  the  multiplicative  congruental1  (power  residue)  method. 

This  method  involves  raultiphylng  two  thirty-five  bit  numbers  together 
(erne  of  them  always  being  the  same)  and  taking  the  least  significant 
thirty-five  bits  of  the  result  as  the  random  number  generated.  This 
number  is  also  used  as  the  next  multiplier.  The  success  of  this 
method  depends  on  the  choice  of  the  constant  and  on  the  choice  of  a 
starting  random  number.  The  method  is  extremely  fast  (about  8 /is  as 
coded  by  the  author  on  a UNIVAC  1108)  and  will  with  proper  normalization 
present  numbers  evenly  distributed  between  0 and  1 with  negligible 
dependence  from  number  to  number. 

The  numbers  are  then  transformed  to  a Gaussian  or  normal  distribu- 
tion with  a zero  mean  and  unit-variance  using  a polyncminal  fit  to  the 
cumulative  distribution  curve.  These  numbers  are  now  equivalent  to 
the  numbers  that  would  be  obtained  after  sampling  a bandlimited  low  pass 
noise  spectrum  of  1 volt  RMS  at  twice  the  cutoff  (Nyquist)  frequency. 

The  remaining  task  is  to  produce  correlation  between  the  numbers  so 
that  they  look  like  numbers  which  were  obtained  from  sampling  a bandpass 
process  at  several  times  Nyquist.  This  is  accomplished  by  sane  form 
of  digital  filtering. 

The  theory  and  application  of  digital  filtering  has  been  documented 
many  places^-*  and  iS  founded  on  Z-transform  theory®.  The  most  efficient 
method  in  terms  of  machine  cycles  is  usually  a filter  of  the  recursive 
type.  Such  an  algorithm  works  with  future,  present  and  past  samples  to 
properly  filter  the  data.  It  was  desired  for  the  initial  simulation 
to  synthesize  a bandpass  filter  in  the  frequency  domain  with  rather 
steep  skirts  so  a 5th  order  Butterworth  design  was  used.  Such  a filter 
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has  the  response  shewn  in  Figure  2 and  needs  17  multiply  adds  per  data 
point.  Using  the  above  techniques  the  author  has  been  able  to  produce 
over  50  channels  of  independent  filtered  noise  about  20,000  samples 
long  in  under  two  minutes  on  the  UNIVAC  1108. 

Synthesis  of  Signals 

This  program  is  designed  to  create  multi-channel  signals  with 
any  desired  characteristics  to  be  added  to  the  synthesized  noise. 
Normally  in  sonar  applications  the  signal  is  distinguished  by  its 
spatial  distribution  being  a plane  wave,  and  by  sane  spectral  shape. 

To  properly  form  the  signal  a time  delay  structure  is  specified  and 
a single  channel  of  filtered  noise,  sinusoidal  components,  etc.  is 
split  and  properly  delayed  into  the  correct  number  of  channels.  Since 
the  signal  is  virtually  formed  on  a sample  by  sample  basis  the  type 
of  signals  formed  may  vary  with  time.  This  allows  simulation  of  fading 
channels,  transients,  and  echoes.  The  time  delay  structure  for  the 
array  under  simulation  and  the  type  of  signal  are  data  inputs  to  this 
program.  The  signal  to  noise  ratios  are  controlled  at  the  time  the 
signal  is  added  to  the  spatially  correlated  noise  produced  in  the  third 
program. 


It  was  desired  to  structure  the  signal  and  noise  for  this  simula- 
tion to  have  a realistic  acoustic  spatial  structure  as  well  as  temporal 
structure.  In  order  to  accomplish  this  a method  for  producing  inter- 
hydrophone  correlation  for  the  noise  field  was  necessary  since  seme 
arrays^*®  are  operating  in  such  an  environment.  The  technique  for 
producing  the  interhydrqphone  correlation  is  detailed  in  a USL  Report 
by  Eby^  and  will  be  reviewed  for  the  reader  in  a simple  form,  as  used 
by  this  program. 

Let  e1(t),  e2(t),  - - - ,eN(t)  represent  the  voltage  at  the  N 

hydrophones  in  question,  having  a specified  covariance  matrix  R = (r^) 
where  rij  = ei(t)ej(t).  The  results  of  the  first  program  produced  N 
independent  channels  of  noise,  which  will  be  represented  as  *i(t),xg(t). 
The  noise  (uncorrelated  between  channels)  covariance  is  Xj.(t)xj (t)  = 0 
for  all  i jt  j.  In  order  to  introduce  correlation  between  the  channels 
an  NXN  matrix  of  real  frequency  independent  weights  A = (aij)»  which 
will  be  called  the  connection  matrix,  must  be  found.  The  specified 
covariance  matrix  may  be  expressed  in  terms  of  the  connection  matrix 


R = AAJ 
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It  turns  out  there  are  many  A matrices,  (at  least  2N)  which  will 
satisfy  (l)  for  a particular  R matrix.  A particular  A matrix  may  be 
found  by  solving  R for  its  N eigenvalues,  denoted  as  x ; , forming  a 
diagonal  matrix  D whose  main  diagonal  consists  of  the  X,  , and  a matrix 
C whose  ith  row  are  the  elements  of  the  ith  eigenvector  associated  with 
R.  Then 

A = D?C  (2) 

where  is  the  square  root  of  the  diagonal  matrix  D. 

Another  method  for  solving  equation  (l)  due  to  Marsaglia10,  provides 
an  algorithm  for  finding  the  elements  of  a triangular  matrix  A in  terms 
of  the  original  covariance  elements  &3  follows: 
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Using  either  of  these  solutions  a network  such  as  the  one  in 
Figure  3,  may  be  used  so  that  a particular  hydrophone  voltage  e^t) 


XI 
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et(t)  = > aJ  . x^(t)  (8) 

is  in  general  a linear  combination  of  the  N uncorrelated  channels. 


Since  the  correlated  noise  field  is  formed  on  a sample  by  sample 
basis  the  interhydrophone  correlation  may  be  programmed  as  a function 
of  time.  In  addition  this  method  allows  the  study  of  arrays  operating 
in  acoustic  noise  fields  that  have  unequal  correlation,  between  equi~ 
distant  hydrophones. 
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The  third  program  is  also  used  to  add  the  signals  generated  in 
the  second  program  to  the  correlated  noise  at  any  desired  s/N  ratio 
at  any  phone.  Multi-target  situations  are  synthesized  by  the  repeated 
addition  of  various  signals,  and  time  varying  S/N  ratios  are  produced 
by  changing  the  S/N  ratios  as  a function  of  time. 

Forming  of  Multiple  Beams 

The  fourth  program  is  designed  to  form  multiple  beams  from  the 
noise  and  signals  that  have  been  synthesized  in  the  first  three 
programs.  A pre-calculated  time  delay  structure  for  the  particular 
array  geometry  of  interest  is  used  to  form  any  number  of  desired 
beams.  In  addition  different  non-linear  operations  may  be  performed 
on  the  data  before  beamforming  so  that  different  systems  may  be  compared 
under  identical  noise  and  signal  characteristics.  The  high  order 
statistics  at  the  beamformer  output  present  considerable  analytical 
difficulties  for  certain  of  the  non-linear  operations,  but  are  of 
considerable  interest  in  system  design  and  performance,  and  may  be 
estimated  from  this  program. 

Post  Beamforming  Processing 

After  beamforming  a variety  of  algorithms  exist  for  multibeam 
sonars  that  deal  with  normalization,  detection,  tracking  and  parameter 
estimation  that  are  easily  programmed.  Simulations  of  these  algorithms 
are  not  meaningful  however  if  the  high  order  statistics  are  unknown  at 
the  inputs,  and  what  their  relationship  is  to  the  input  acoustic  field. 
Therefore  they  must  be  an  integral  part  of  the  total  simulation.  An 
example  of  the  difficulties  encountered  is  shown  by  considering  the 
problem  of  trying  to  locate  an  echo  on  three  adjacent  beams.  What  is 
the  probability  of  an  echo  being  detected  on  two  of  the  three  adjacent 
beams?  Do  false  alarms  occur  at  the  same  time  on  all  three  beams? 

The  problems  are  clearly  of  a multidimensional  nature,  where  the 
efficiency  of  Monte  Carlo  techniques  becomes  competitive  with 
theoretical  methods. 

Statistical  Estimation  and  Utility  Routines 

The  last  main  program  is  designed  to  provide  estimates  of  the 
relevant  statistics  at  various  points  in  the  algorithm  or  "system". 
Various  plots  and  graphs  are  produced  and  summaries  of  each  run.  In 
addition,  various  routines  are  used  for  tape  handling,  since  tapes 
are  produced  between  each  program  so  that  different  algorithms  may  be 
compared  with  the  same  inputs.  Sane  of  these  routines  have  been 
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detailed  in  other  USL  Technical  Memorandum 


Conclusions 


Monte  Carlo  simulations  can  provide  detailed  knowledge  of  how 
major  and  minor  changes  in  an  algorithm  (and  the  corresponding  special 
purpose  hardware  used  to  implement  the  algorithm)  effect  certain  sonar 
performance  parameters,  under  Identical  and  controllable  input  conditions. 
Because  of  this  it  provides  a unique  sonar  design  capability  that  may 
be  effectively  used  in  trade-off  studies.  In  most  cases  the  simulations 
run  on  a general  purpose  computer  will  not  run  as  fast  as  the  special 
purpose  hardware,  and  therefore  should  not  be  used  in  a major  data 
reduction  program. 
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Figure 


Input  Noise  Spectrum  as  Synthesized  by 
Digital  Filtering  of  Itocorrelated  Noise  Samples 


